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Abstract 

Our quasi-particle model is compared with recent lattice QCD data at finite 
temperature and baryon number density with emphasis on the coefficients in 
q ■ the Taylor series expansion of thermodynamic observables. The inclusion of 

CN) . static critical end point effects into the equation of state is discussed. 
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1 Introduction 

The QCD phase diagram exhibits an astonishingly rich phase structure. The (pseudo-) critical line, which 
separates the phase dominated by quark and gluon degrees of freedom at large temperature T and baryo- 
chemical potential /xg from the one dominated by hadrons and resonances, has been investigated by 
means of lattice gauge theory Q|2l|3]|4]]. From theoretical reasoning [5 6 1 for finite quark masses, the 
phase transition is of first order at finite T and large /j, b ending in a critical point of second order. For 
smaller thermodynamic observables change rapidly but continuously indicating a crossover regime. 
There, the equation of state (EoS) has been computed for Nj = 2 quark flavours QUI- While the location 
of the critical end point (CEP) showing a strong quark mass dependence |3j was determined by first 
principle QCD evaluations 000), the extension of the critical region is fairly unknown. Lattice QCD 
studies of the volume dependence of the Binder cumulant indicate that the CEP belongs to the 
static universality class of the 3D Ising model. At present, many investigations aim to study implications 
of such a fundamental issue of QCD. In particular, observable consequences of the occurrence of the CEP 
as novel feature of QCD are discussed [ 12 1. Heaving in mind the successful hydrodynamical description 
of the expansion stage in heavy-ion collisions, one intriguing problem concerns the manner the EoS of 
strongly interacting matter becomes modified by the CEP. This is the subject of the present contribution. 

In |2j the quasi-particle model is shortly reviewed and compared with recent lattice QCD results. 
We focus on the Taylor series expansion coefficients QUI. In|3j CEP effects on the EoS are discussed 
for a toy model and for our QCD based quasi-particle model. The results are summarized in|4j 

2 Taylor series expansion of the EoS 

Thermodynamic observables can be expressed as Taylor series expansions in powers of /j,b/T. Accord- 
ingly, the pressure is decomposed into 

p(T,^) = r 4 ^ C2 „(r)^j , CD 

where c (T) = p(T,n = 0)/T 4 and c k {T) = d k p/dp k T k ~ 4 /kl with n = fi B /3. In jTHH, the 

Taylor expansion coefficients up to cq(T) have been presented basing on first principle QCD evaluations. 

Achieving a flexible parametrization of the EoS, we formulated a model which describes the quark- 
gluon fluid in terms of quasi-particle excitations fT3l[T4l 

p(T,fj)= Pi(T,p,)-B(T,p,). (2) 
i = g,9 
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Fig. 1: Taylor expansion coefficients co,2,4 (top 
row, horizontal lines indicate corresponding co- 
efficients for a non-interacting gas above T c 
as previously often used to describe the quark- 
gluon plasma) and cs (left panel in bottom row) 
of as a function of T/T c . Lattice QCD data 
(symbols) from (7||8|. 
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Here, p. L denote thermodynamic standard expressions for quarks and transverse gluons with dynamically 
generated self-energies Lfj and a non-perturbative effective coupling G 2 (T, p) as essential input. Ther- 
modynamic self-consistency is ensured through the stationarity conditions Sp/SILi = 0, imposing in turn 
conditions onto B(T, p). From Maxwell's relation for p, a flow equation for G 2 (T, p) follows [ 13 1 



dG 2 dG 2 
+ ax- 



dp 



dT 



b. 



(3) 



Knowing G 2 on an arbitrary curve in the T - p plane, (J3} can be solved as a Cauchy problem. For 
convenience, we adjust G 2 (T, p = 0) to lattice data at p = enabling a mapping into the finite chemical 
potential region via ©. The Taylor expansion coefficients follow straightforwardly from © as integral 
expressions involving G 2 and higher order derivatives of the effective coupling at vanishing chemical 
potential. The latter can be computed by exploiting the flow equation © (cf. [ 14] for details). In Fig. 1, 
a fairly good agreement between Taylor expansion coefficients from the quasi-particle model and lattice 
results is shown. In particular, the peak in c±(T) and the dipole structure in cq(T) are reproduced. 
Adjusting the parametrization of the effective coupling onto cq{T) dictates a change in G 2 (T,p = 0) 
at T c = 170 MeV from a regularized logarithmic dependence (resembling the perturbative behaviour at 
large T) into a linear dependence. This change in the curvature of G 2 can be considered as implemented 
phase transition and is responsible for the pronounced structures in c±$(T). In fact, these structures 
disappear when neglecting higher order derivatives of the effective coupling in the integral expressions 
°f Ci,e{T) serving for a test of the flow equation (J3J- 



3 Critical end point 

Stalling from a thermodynamic potential, e. g. Gibbs free enthalpy, it can be decomposed into an analytic 
and a non-analytic part where the latter is related to phase transitions and critical phenomena 1151 . Ac- 
cordingly, the EoS formulated in terms of the entropy density is given through s = s a + s n . Here, the an- 
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Fig. 2: Isentropic trajectories of the toy model J4j depending on the strength parameter A for s/ns = 50, 28. Dashed, thin 
and solid lines exhibit results for A — 0., 0.5, 1.0 respectively. Dotted lines represent the tangent on the estimated (pseudo-) 
critical line at the CEP. Left panel: AT = 100 MeV, A^ s = 200 MeV, D = 0.15 ; right panel: AT = 10 MeV, A^ s = 10 
MeV, D = 0.06. 



alytic contribution s a has to be adjusted to the known EoS outside of the critical region. The non-analytic 
part s n should embody the feature of being continuous left to the CEP (i. e. at small chemical potentials) 
whereas on the right (at large chemical potentials) it generates a first order phase transition. A convenient 
parametrization of s n for the 3D Ising model characterized by a set of critical exponents is given in [ 16 1. 
Still, the corresponding variables employed usually in condensed matter physics including the order pa- 
rameter need to be mapped into the T - \ib plane in the vicinity of the CEP. Details of this mapping and a 
useful formulation of the entropy density contribution can be found in the pioneering work FPTll we rely 
on. In the following, we estimate the phase border line to be given by T c (/xb) = T c (l + ^d(/is/3T c ) 2 ^ 
with d = —0.122 according to [ 1 3 1 and locate the CEP at \i B ,c = 360 MeV in agreement with J4). 

As a simple toy model, let us employ the first terms in Q, however, with constant expansion 
coefficients. The entropy density contributions are given by 

2 2 
s a (T,n B ) = 4c T 3 + -c 2 /i|T, s n (T,fi B ) = -c 2 n 2 B T A tanh (S c (T,fi B )) (4) 
y y 

with c = (32 + 21A r / )7T 2 /180, c 2 = Nj/2 and N f = 2. n B follows from © via standard thermo- 
dynamic relations (cf. fTTlO with integration constant ne(0, /j,b) = §C4(//b/3) 3 where C4 = Nj/A-k 2 . 
The ansatz for s n has been chosen such that s n — > for T — > and the net baryon density vanishes at 
[is = 0. The parameter A describes the strength of the non-analytic contribution in the EoS. We apply 
the same S c {T,hb) as in lUTl assuming a fairly large critical region parametrized by AT = 100 MeV, 
A/ib = 200 MeV and a stretch factor D = 0.15. Hence, CEP effects on the EoS and in particular on 
isentropic trajectories s/ns = const in the T - (j,b plane can be demonstrated. In Fig. 2, the influence of 
the strength parameter A on the behaviour of isentropic trajectories is exhibited. For increasing A > 0, 
the trajectories for large s/tib tend to be attracted towards larger u b due to the presence of the CEP. In 
fact, the CEP acts as an attractor on trajectories on the left whereas on the right a repulsive impact is 
found. Evidently, the curves on the right side of the CEP display the existence of the first order phase 
transition. By shrinking the critical region to AT = 10 MeV, Au B = 10 MeV and D = 0.06, the 
influence of the CEP decreases (right panel of Fig. 2) in comparison with the results obtained for a large 
critical region (left panel of Fig. 2). In particular, the sections on the hadronic side become less affected 
when decreasing the extension of the critical region. 

The parameters in the non-analytic entropy density contribution and in particular A have to be 
chosen such that standard thermodynamic consistency conditions are satisfied iTPTl . Accordingly, during 
the adiabatic expansion of the system and its related cooling, both n b and s must decrease. For A < 
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Fig. 3: Isentropic trajectories in the truncated quasi-particle model adjusted to lattice QCD |7][8) without CEP (left panel) 
and with CEP effects included (right panel) parametrized by AT = 10 MeV, Afi B = 10 MeV, D = 0.06 and A = 0.5 for 
s/ub = 200, 100, 50, 33 (from left to right). Dotted line represents the estimated phase border line, long-dashed lines in the 
left panel depict resonance gas trajectories with corresponding chemical freeze-out points from 1 18 1. 

with trajectories for large s/ng bent to smaller hb due to the CEP inclusion, however, these conditions 
are violated in the vicinity of the first order transition line. Therefore, the pattern of the isentropic 
trajectories as exhibited in Fig. 2 seems to be generic (cf. fTTt for a different model with CEP which 
involves a hadronic low-temperature and a partonic high-temperature phase). Clearly, this statement 
decisively depends on s in the hadronic phase where the simple toy model cannot account for QCD. 

In contrast, being inspired by lattice QCD we construct the EoS as truncated Taylor series ex- 
pansion including the coefficients cq,2,4,g(T) of Fig. 1. These lattice based trajectories are shown in 
Fig.|3](left panel) where the pattern differs notably from the observations made in the above toy model. 
Although first principle evaluations still suffer from too heavy quark masses and deviate consequently 
from the resonance gas trajectories in the low-temperature phase, the values s/ng on these lattice QCD 
deduced curves agree with the values at the according nearby chemical freeze-out points inferred from 
data [ 18 1. Furthermore, in the deconfined phase lattice results are unstable due to convergence radius 
studies [8 1 at least up to \i B = 300 MeV. Therefore, adjusting the analytic contribution of the EoS known 
from the lattice data by our quasi-particle model, critical end point effects should become visible only for 
larger \i b implying a small critical region. We include the CEP in line with the procedure outlined above 
in our quasi-particle model replacing C2 by C2 (T) in s n and considering a small critical region character- 
ized by AT = 10 MeV, A/j, B = 10 MeV and D = 0.06 with strength parameter A = 0.5. As exhibited 
in Fig. |3] (right panel), CEP effects on isentropic trajectories are significant only for large /j,b with neg- 
ligible impact on the hadronic sections. Nonetheless, the baryon number susceptibility \B = d 2 p/dfi B 
being a measure for baryon number fluctuations diverges for \i B > Hb,c (Fig. 4 right panel) due to the 
discontinuity evolving in ub in contrast to the analytic behaviour (Fig. 4 left panel) stemming from the 
quasi-particle model not containing CEP effects. 

4 Conclusion 

Our quasi-particle model without implemented CEP was successfully compared with recent lattice QCD 
results of the Taylor series expansion coefficients cq(T) and C2,4fi(T). Accordingly, our extrapolation 
procedure into the finite chemical potential region was tested. We considered simple models including 
phenomenologically the QCD critical end point and studied the effects on isentropic trajectories. We 
followed [17] and looked for indications of the CEP acting generically as attractor or repulsor. In fact, 
this is of interest with respect to the question whether CEP effects show up in heavy ion experiments only 
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Fig. 4: Scaled baryon number susceptibility of the quasi-particle model neglecting C6 (T) without CEP (left panel) and with 
CEP inclusion (right panel) as function of T/T c for \i B = 450, 330, 300, 150 MeV (from top to bottom). 

in a very narrow beam energy range. Clearly, appropriate dynamical simulations are needed to account 
properly for such questions. The study of the pattern of isentropic trajectories is only a first step towards 
elucidating possible implications of the very existence of the CEP in QCD. 

Inspiring discussions with M. Asakawa and F. Karsch are greatfully acknowledged. The work is 
supported by BMBF, GSI, EU-I3HP. 
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